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Abstract 

A spherically symmetric fluid code that includes the relativistic efl'ects of hot electrons has 
been newly developed. In the present simulations, I attempt to implode a target shell of laser 
fusion using a nanosecond pulse; then, an additional heating short pulse is injected into the 
coronal plasma surrounding a highly compressed core. I found that the relativistic effects of hot 
electrons on electron transport are likely to inhibit the heat flows, and to reduce thermonuclear 
fusion neutrons significantly. This feature may be important for off-center fast ignition and burn 
of fusion targets. 
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The relativistic effects of liot electrons in laboratory plasmas have attracted much inter- 
ests for the past decades , particularly, relating to laser fusion fast ignition (LFFI) whose 
scenario is achieved by standard high-compression, irradiation of a hole-boring laser beam, 
irradiation of an ignitor laser beam, and thermonuclear burn P| . When high- intensity short 
pulse laser light is impinged on the coronal plasma surrounding a highly compressed core in 
the 'hole-boring' and 'ignition' stages, energetic particles are sprayed through various mech- 
anisms. In order to investigate the relativistic electron transport carrying huge currents, 
substantial studies were performed on the multidimensional relativistic particle-in-cell as 
well as magneto-hydrodynamic jj] simulations. 

Moreover, the fast ignition conditions of the compressed targets have been presented 
in some publications jJO, Q, Ol- Indeed, it is important to survey the energy transport 
and relaxation in dense plasmas over a large spatiotemporal system of R > 10^ — 10^ fim 
and t > 1 ps. Concerning the extension of the analysis of ignition conditions, we are now 
strongly encouraged to irnprove the transport coefficients so as to compare them with those 
for kinetic simulations 0,1^ and experiments H, 3|- That is why, I attempted to construct 



a computational fluid code that includes relativistic effects of electron transport. In this 
code, the transport properties of hot electrons are fully consistent with the current-neutral 
electric flelds self-induced in plasmas, as discussed in ref. |[l3j in detail. 

In this paper, I argue that in the context of the LFFI, heat flux inhibition owing to 
relativistic effects of electrons degrades the transport efficiency of deflagration thermal waves, 
and leads to some reductions in the neutron yields. This is essentially because the drift 
velocity carrying heat asymptotically approaches the speed of light. Since the velocity 
moment of the heat flux is large, the energy of electrons carrying dominant heat (EE-CDH) 
is much larger than the thermal energy: Tee-cdh ~ 7Tth. Thus, it is expected that the 
relativistic effects of electron heat transport can be observed even in the energy domain 
of Tth ~ 10^ — 10^ keV. It is noteworthy that such situations can be easily established 
by irradiating a moderate- intensity laser light of /lA| > 10^^ W/cm^ — /im^ \14\. As for 
astrophysical aspects, the issue presented here is relevant to supernova outbursts l^ . 

The spherically symmetric fluid code has been developed on the Lagrangian frame dm = 
Anpr'^dr with the speciflc volume of = p~^. This solves a set of two-temperature fluid 
equations for an ideal gas, including the relativistically corrected equations for the internal 
energy of electrons. That is. 
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Here, conventional abbreviations have been used. In particular, K,e.rei, ^ei.rei, and c^e^rei 
denote the relativistic Spitzer-Harm heat conductivity energy relaxation coefficient 



jlpj l , and specific heat 
values are given by 



17| . respectively. These coefficients normalized by their nonrelativistic 
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where Ky{a) is the modified Bessel function of index v with its argument of a = mQC^/Tg, 
and the functions 0i(a), 62(a), and 63(a) are defined by [3] 
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where Ei(— a) is the exponential integral function. The normalized coefficients (^J-® for 
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some a values are shown in Table HI Note that for a 1, eq. (jlj) asymptotically approaches 
unity, and for a <^ 1, approaches the expression of Re = [5(27r)^/^/32]a;^/^ |l8|. Another 
important point is that the heat capacity tends to increase in high-temperature regimes, up 
to twofold in the ultrarelativistic limit. Making use of the implicit finite differential scheme, 
we can numerically integrate the internal energy eqs. ©-Q as outlined in the Appendix. 

In the LFFI context, of course, one can also construct some kinds of kinetic simulation 
code. Q to inve.,gate a.y..et* .an.po. ,.voM„, magnetic fields MQ. but the. 
allowed spatiotemporal ranges are small. Regarding the present scheme, the macroscopic 
transport properties in supersolid density regions and the effects of a more realistic plasma 
gradient as well, can be revealed as a trade-off for missing the smaller spatiotemporal scales, 
such as the Debye length and the plasma oscillation period. Along the diffusion approx- 
imation adopted under geometrical constraint, irradiating laser intensity II is traced by 
v.iaiJOr) ^ where^and ... denote the group velocity of hght and the absorp- 

tion coefficient, respectively [20|. For simplicity, in the present simulations, the laser light 
is set at normal incidence. The ray deposits its own energy, when propagating through the 
corona, and is then resonantly damped at the critical density. Note that in eq. (0) radiation 
pressure Pph is also taken into account. 

Furthermore, we are concerned with the thermonuclear reactions of deuterium-tritium 
(DT) and deuterium-deuterium (DD) fusion. The rate equations for number density of the 
tritium Nt and deuterium Nd |2L] were solved in a postprocessing manner. Assuming the 
local self-heating due to charged particles, one obtains the source terms of eqs. (j2I) and (jH)) 
as functions of the heating power densities of Sdt oc {vaoT) N^Nt and Sdd oc {vgod) 
[2^ . Although these models are rather crude, they still have merit for the investigation of 
fundamental transport processes of electrons. As shown later, it is instructive to compare 
the neutron yields derived from various models of electron transport. 

In the following, I demonstrate how a highly compressed fuel can be heated due to 
electron transport, and ignited. For convenience, 'case A' and 'case B' are used as references 
for simulations with and without relativistic effects, respectively. The initial laser conditions 
and target parameters are shown in Table |n] Hereafter, the target parameters are fixed. 
The total mass of the target shell is on the order of that presented in ref. |23], and the initial 
aspect ratio of the shell is about {R/AR)q ~ 10. High compression of the shell, maintaining 
a low entropy, is carried out by the Gaussian pulse shaping of driver laser light. The outer 
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thin ablator of carbonized DT (CDT) is blown off just before the deceleration phase, leaving 
the dense compressed core of DT fuel. Hydrodynamic instabilities concomitant with the 
low-entropy implosion are omitted for the moment. 

In Fig. n^a) for case A, I show a flow diagram of implosion with additional heating. 
The total mesh number is J = 102, and the flow lines are displayed each tenth mesh. 
At t = 2.25 ns, an additional 10 PW (10^^ W) power with the pulse duration of r = 
1 ps is deposited at a relativistically modified cut-off density which is defined by nc_rei = 
ric^/l + (1.37 X 10^8 W/cm2 - fxra^), where ric = 9.97 x lO^o (1.06 iim/Xif cm'^ is 

the nonrelativistic cut-off density. In this case, the pulse intensity becomes about /lA| ~ 
10^8 W/cm^ — /im^. One can see of the detonation shock propagating radially inwards, as 
well as the explosion, associated with the off-center fast ignition. The flow diagram without 
additional heating is shown in Fig. ^b). It is confirmed that ignition does not occur in the 
standard implosion of the small target when using a driver laser energy of about 4 kJ. 

The spatial profiles of plasma temperature and mass density at t = 2.251 ns, the moment 
that the additional heating power of 10 PW is switched off, are shown in Fig. |21for cases A 
and B. One may notice that in the regions of r = — 35 yum, the stagnating isobaric 
core where a central hot spot is surrounded with cold dense plasma (DT main fuel) is well 
established as a result of the low-entropy implosion. The equilibrium static pressure reaches 
P = Pe + Pj > 10 Gbar. It turns out that nonlinear propagation of the deflagration wave 
directly heats the main fuel. Here we define the average electron temperature of the main 
fuel as (Te) = J^^ pT^r'^dr / J^^ pr'^dr, where ri and r2 indicate the inner and outer radii 
where p = Ppeak/10, respectively. In case B without relativistic corrections, the average 
temperature is (Tg) ^ 9.9 keV. On the other hand, in case A, the heat flux inhibition due 
to relativistic effects lowers the temperature to (Tg) ~ 6.6 keV. This reflects the factor 
ke/cve — 0.75 for Tg ^ 25 keV in the tenuous corona (see Table |ll). The relativistic effects 
on energy transfer between electrons and ions are small. 

In Fig. 121 for case A, I show the spatial profiles of plasma temperature and mass density, 
just before the irradiation of the ignitor pulse {t = 2.25 ns), during the irradiation [t = 
2.2505 ns), and just after the irradiation {t = 2.251 ns). While the coronal electrons are 
rapidly heated up to Tg ~ 25 keV, the deflagration wavefront slowly propagates in the cold 
dense region where the heat capacity is very large, and the steep temperature gradient at the 
wavefront self-generates a longitudinal electric field having a maximum value on the order 
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of magnitude of i?max ~ e~^\dT(./dr\ ~ 10^" V/m [18]. For the plasma parameters shown in 
Figs. 121 and 01 the heating rate of the DT fuel ions is estimated to be dTj/dt ~ 1 — 1.5 keV/ps, 
when assuming the Coulomb logarithm of InA ~ 5. The fuel can be, therefore, heated to 
the ignition temperature until the stagnating core is disassembled. 

Let us show the additional heating power dependencies of neutron yields in Fig. |^ for 
both cases A and B. Within the range of 0.1 — 20 PW, an additional power with the pulse 
duration of r = 1 ps is injected at t = 2.250 ns. It is found that in the present model the 
threshold power for ignition seems to be about 1 PW, corresponding to the energy of 1 kJ. 
The pulse power of about 20 PW leads to a considerable increase in neutron yields, which 
are of the order of one thousand times that in the case without additional heating. The 
time-integrated bremsstrahlung loss is 60 — 70 J at most. It is noted that the flux inhibition 
degrades the heating efficiency, thereby reducing the neutron yields, particularly, when the 
additional power is, in this model, in the range of 3 — 10 PW. The yields tend to be sensitive 
to the timing of the irradiation of the ignitor pulse 0, |3| • 

The present predictions of the neutron yields and the threshold power for ignition may be 
rather pessimistic for determining the lower limit of the yields. In order to argue this point, 
I show, in Fig. |2 some results of the Fokker-Planck simulation coupled with spherically 
imploding hydrodynamics In this calculated example, the laser irradiation conditions 
and target parameters are similar to those presented in Table |n] It is assumed that the 
absorbed energy of an ignitor pulse is transfered by 50 % into the high-energy tail electrons 
of the temperature (variance) of 50 keV. As shown in the figure, the effects of nonlocal 
transport seem to increase the neutron yields, since the hot tail electrons, which possess a 
longer mean-free path, deeply penetrate into the compressed plasma, preheating the fuel. 
It is noted that such nonlocality, as well as anisotropy of energy deposition, could be an 
advantage to the LFFI. The details should be clarified by multidimensional fiuid simulations 
including relativistic kinetics in the future. 

APPENDIX: IMPLICIT DIFFERENCE OF RELATIVISTICALLY COR- 
RECTED INTERNAL ENERGY EQUATIONS 

In this appendix, I briefiy explain the numerical method to integrate the relativistically 
corrected energy equations. The implicit finite difference of eq. ^ can be cast to 
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— (1 — 0) (x'gj ,,-_^i/2 (^e,i+l/2 ~ ^"-+1/2) ~ 0^ri,j+l/2 (^e"j+l/2 ~ Tl^j+l/2) + '^'"+1/2; (^-1) 

where (j = 0, 1, 2, J) denotes the discrete positions of fluid elements and (< 1/2) is 
the implicit parameter, e.g., the Crank-Nicholson scheme corresponds to the case of = 1/2. 
The last term on the right-hand side, 5'"^^^25 represents generic power sources and/or sinks. 
Note that the coefficients He, ^ei, and Cye depend upon the temperatures. 

In the present simulation, the tridiagonal matrix (A.l) for = is inverted by the cyclic 
reduction method [2^. The temperature at the intermediate time steps T*j_^_^^^ and the 
coefficients involving this are iteratively advanced to the next {T**_^i^2) ■ When the condition 
of |T;*.^i/2 - ^e,j+i/2l/ (^ej+1/2 + ^ej+1/2) < e (y j) is satisfied for a small value of e, the 
vector Tg*^j^y2 is replaced by These procedures are repeated each time step, whose 

increment At"+^/^ is automatically changed, invoking the Courant-Friedrich-Lewy condition. 
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(3) Radius r(^im) (b) 

FIG. 1: Flow diagrams of laser implosion (a) with additional heating and (b) without additional 
heating. The figures have the same axes, and the flow lines are displayed each tenth mesh. At 
t = 2.250 ns for (a), the additional heating power of 10 PW having the pulse duration of r = 1 ps 
is deposited at the relativistically corrected critical surface. 
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FIG. 2: The spatial profiles of electron temperature (full curves), ion temperature (dot-dashed 
curves), and mass density (dotted curves) at t = 2.251 ns (a) for case A with relativistic effects 
and (b) for case B without relativistic effects. The figures have the same axes. 
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FIG. 3: The spatial profiles of electron temperature (full curves), ion temperature (dot-dashed 
curves), and mass density (dotted curves) for case A; (a) just before the pulse power injection, (b) 
during the injection of 10 PW, and (c) just after that. The figures have the same axes. 
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10"" 10° 10^ 10^ 

Additional heating power (PW) 

FIG. 4: Additional heating power dependencies of neutron yields. The full and open circles are for 
cases A and B, respectively. The target parameters and driver laser conditions for implosion are the 
same as those given in Table ^1 For comparison, the corresponding results of the Fokker-Planck 
(FP) simulation are also shown (open squares). For explanation see text. 
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TABLE I: The relativistic correction factors given by eqs. ©-(jTl) in the text. 



a 


Te (MeV) 






'-ve 




< 1 


— 


— 


— 


2 


— 


0.05 


10.2 


0.087556 


5.6084 


1.9992 


0.043796 


0.1 


5.11 


0.12373 


3.9726 


1.9968 


0.061966 


0.5 


1.02 


0.27123 


1.8510 


1.9393 


0.13986 


1 


0.511 


0.36792 


1.4188 


1.8343 


0.20058 


5 


0.102 


0.62668 


1.0528 


1.3746 


0.45589 


10 


0.0511 


0.73598 


1.0206 


1.2156 


0.60546 


20 


0.0256 


0.83044 


1.0084 


1.1160 


0.74420 


» 1 




1 


1 


1 


1 



'^Cited from ref. 
^ Cited from ref. 
Cited from ref. 
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111 



by Honda and Mima. 

by Beliaev and Budlcer. 

by Balescu and Paiva-Veretcnnicoff. 
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TABLE II: Simulation parameters of laser implosion including additional heating. 





Laser conditions 
Driver pulse 


Ignitor pulse 


Target parameters 


Pulse shape: 


Gaussian 


Square 


Inner radius: 


223.32 /im (vacuum) 


Laser energy: 


4.1 kJ 


10 kJ'' 


D0.5T0.5: 


2.00 ^m (3 meshes) 


Wavelength: 


0.53 ^m 


1.06 ^m 


D0.5T0.5: 


16.92 fim (42 meshes) 


Rise time: 


0.6 ns* 




C0.426D0.534T0.0153 : 


4.76 ^m (57 meshes) 


Pulse width: 


1.909 ns'= 


1.0 ps'^ 






Fall time: 


0.5 ns*^ 








Peak power: 


2.5 TW 


10 FW 







''Variable between 0.1 kJ-20 kJ (100 TW-20 PW), as seen in Fig.H 
'it follows the Gaussian pedestal of 1.2 ns. 
'^FWHM, including the top flat part of 0.809 ns. 
''Switched on at i = 2.250 ns. 

•^Additional Gaussian pedestal of 1.0 ns follows its fall. 
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